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Abstract Evidence shows that massive black holes reside in most local galaxies. Stud- 
ies have also established a number of relations between the MBH mass and properties 
of the host galaxy such as bulge mass and velocity dispersion. These results suggest 
that central MBHs, while much less massive than the host (~ 0.1%), are linked to the 
evolution of galactic structure. In hierarchical cosmologies, a single big galaxy today 
can be traced back to the stage when it was split up in hundreds of smaller compo- 
nents. Did MBH seeds form with the same efficiency in small proto-galaxies, or did 
their formation had to await the buildup of substantial galaxies with deeper potential 
wells? I briefly review here some of the physical processes that are conducive to the 
evolution of the massive black hole population. I will discuss black hole formation pro- 
cesses for 'seed' black holes that are likely to place at early cosmic epochs, and possible 
observational tests of these scenarios. 
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1 Introduction 



Black holes, as physical entities, span the full range of masses, from tiny holes pre- 
dicted by string theory, to monsters as massive as a dwarf galaxy. Notwithstanding 
the several orders of magnitude difference between the smallest and the largest black 
hole known, all of them can be described by only three parameters: mass, spin and 
charge. Astrophysical black holes are even simpler systems, as charge can be neglected. 
Complexities arise because of the interaction between astrophysical black holes and 
their environment. 
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I will focus here on the formation and evolution of massive black holes (MBHs), in 
high-redshift galaxies, and their symbiotic evolution with their hosts. MBHs weighing 
million solar masses and above have been recognized as the engines that power quasars 
detected at early cosmic times. Dynamical evidence also indicates that MBHs with 
masses in the range M^h ~ 10^ — 10^ ordinarily dwell in the centers of most nearby 
galaxies ( [Ferrarese and Ford[|2005| . The evidence is particula rly compelling for our own 



2003 



Ghez 



galaxy, hosting a central MBH with mass ~ 4 x 10° Mq (e.g., ISchodel et al 
[et al,^2005J . MBHs populate galaxy centers today, and shone as quasars in the past; the 
quiescent MBHs that we detect now in nearby bulges would be the dormant remnants of 
this fiery past. Dynamical estimates indicate that, across a wide range, the central black 



hole mass is about 0.1% of the spheroidal component of the host galaxy (Magorrian 



et al 1998 Marconi and Hunt 2003 Haring and Rix 2004 1 . A tight correlation is also 



observed between the MBH mass and the stellar velocity dispersion of the hot stellar 
component ([Ferrarese and Merritt} |2000| [Gebhardt et "al} |2000| [Tremaine et"ai| |2002| 



Giiltekin et al 20091. The surprisingly clear correlations between MBH masses and 



the properties of their host galaxies suggest a single mechanism for assembling MBHs 
and forming galaxies. The evidence therefore favors a common root, a co-evolution, 
between galaxies and MBHs. These correlations may well extend down to the smallest 
masses. For example, the dwarf Seyfert 1 ga laxy POX 52 is thought to contain a MBH 
of mass Mbh 



W^Mq 



( Barth et al 2004 1 



At the other end, however, some powerful quasars have already been detected at 
2 > 6, corresponding to a time less than a tenth of the age of the Universe, roughly 
one billion years after the Big Bang. Follow-up observations confirmed that at least 
some of these quasars are powered by super- massive black holes with masses ~ 10^ M© 
(Barth et al 2003 Willott et al 2005 I, probably residing in the centers of substantial 
galaxies. However, these exceptionally bright quasars are most likely just the tip of the 
iceberg: rare objects- on the tail of the mass distribution. This implies larger numbers 
of less exceptional objects, and that MBHs existed in large numbers during the Dark 
Ages, before or around the time when the first stars formed. 



We are therefore left with the task of explaining the presence of MBHs when the 
Universe is less than 1 Gyr old, and of much smaller MBHs lurking in 13 Gyr old 
galaxies. The outstanding questions concern how and when "seed" MBHs formed, the 
frequency of MBHs in galaxies, and how efficiently MBH seeds grew in mass during the 
first few billion years of their lives. The formation of MBHs is indeed far less understood 
than that of their light, stellar mass, counterparts, end-points of stellar evolution for 
stars more massive than ~ 20 M© . The "flow chart" presented by Rees ( 1978 1 still 
stands as a guideline for the possible paths leading to formation of massive MBH seeds 
in the center of galactic structures. In the following I will review the main physical 
processes thought to influence MBH formation. I will mostly focus on astrophysical 
processes that happen in galaxies, and I will consider three possibilities: (i) that MBHs 
are the remnants of the first generation (PopHI) stars, (ii) that MBHs formation is 
triggered by gas-dynamical instabilities, (iii) that MBH seeds are formed via stellar- 
dynamical processes. A simplified scheme that describes the possible routes to MBH 
formation in high-redshift galaxies is shown in Figure [l] I will also briefly mention the 
possibility that MBHs are related to inflationary primordial black holes. 
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Fig. 1 Scheme of possible pathways to MBH formation in high-redshift galaxies. Artwork by 
B. Devecchi. 



2 Massive black hole formation 



2.1 Cosmological background 



The observation of luminous quasars at 2: ~ 6 (e.g, Fan et al 2001b I implies that the 



first MBHs must have formed at very early times. The luminosities of these quasars, 
well in excess of 10^^ erg s~^, imply MBHs with masses ~ 1O^M0 already in place 
when the Universe is only 1 Gyr old. The accretion of mass at the Eddington rate 
causes a black hole mass to increase in time as 



M{t) = M(0) exp 



1 



t 



^Edd 



(1) 



where t^dd ~ 0.45 Gyr and e is the radiative efficiency. For a 'standard' radiative 
efficiency e ~ 0.1, and a seed mass M(0) — 10^ — IO^Mq, it takes at least 0.5 Gyr to 
grow to ~ lO^Af0. This brings us back in time to an epoch when the Universe was 
very young and galaxies in their infancy. 

Current theories of structure formation favour their growth from gravitational col- 
lapse of small perturbations in a quasi-homogeneus Universe, dominated by Cold Dark 
Matter (CDM). In the framework of CDM models the collapse of structures proceeds 
bottom-up on larger and larger scales, giving rise to a hierarchy of smaller structures 



that are incorporated into larger ones at later times (White and Rees 1978 1. In the fol 



lowing I briefly summarize the basic points necessary to describe the cosmic landscape 
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where MBH formation takes place. For a detailed description of the cosmic structure 
formation physics, I refer the reader to Peacock (19991. 

The evolution of perturbations up to recombination and on large scales, much larger 
than the typical correlation length of the matter field, is relatively easy to follow. Such 
perturbations remain small all the time, and it is a reasonably good approximation 
to describe them as if they were in linear regime. After recombination, fluctuations 
continue growing, until they turn non-linear. At this stage structure formation takes 
place: collapsed clumps form, which eventually reach an equilibrium configuration at 
virialization. 

By means of the virial theorem and conservation of energy the properties of virial- 
ized halos can be expressed as a function of the virial overdensity with respect to the 
background density, iSvir ~ 178 in an Einstein-de Sitter Universe. Detailed calculations 
for different cosmologies (e.g., /ICDM) can be found in iLacey and Cole (19931 and 



Eke et al (19961. A halo at redshift z is therefore uniquely characterized by a virial 



radius rvir, defined as the radius of the sphere encompassing a mean mass overden- 
sity (5vir. From the virial theorem the virial mass, Mvir can be calculated straightfor- 
wardly, along with the circular velocity, Vc ~ \/ GM^\^lr^\^ and virial temperature 
Tvir = ^.TUpVc / {"ikB) , where /i is the mean molecular weight, mp is the proton mass 
and kg is the Boltzmann constant. 

It is within these dark matter halos growing out of small primordial density fluctu- 
ations by gravitational instabilities that galaxies form out of baryons. The cosmological 
Jeans mass, Mj, represents the minimum mass of an overdense region where pressure 
gradients are sufficiently weak and the gravitational collapse of the baryon component 
can proceed. The first bound systems with M > Mj capture baryons and during sub- 
sequent collapse and virialization of the dark matter the baryons are shock heated to 
the virial temperature. The requirement for a dark matter -I- baryon system to become 
a galaxy is that star formation takes places. Star formation can occur if cooling allows 
the baryonic clouds to dissipate their kinetic energy, continue to collapse and fragment. 

For low-mass objects, the smaller they are, the less efficiently they dissipate en- 
ergy and cool. Thus a detailed treatment of gas-dynamical processes will predict a 
characteristic mass scale such that more massive objects can cool rapidly, whereas 
smaller lumps will merely remain pressure-supported and do not form luminous ob- 
jects (Tegmark et al 19971. The first collapsing halos have virial temperatures smaller 
than the 10 K temperature at which cooling from electronic excitation of atomic hy- 
drogen becomes effective. In order for their gas to cool down and form the first stars, 
these halos must rely on the less effective H2 cooling. Halos with T^ir < 10^ K cooling 
by H2 are typically referred to as 'minihalos'. 



2.2 Population III remnants 

One of the most popular scenarios for MBH formation associates MBH seeds with the 
remnants of the first generation of stars, formed out of zero metallicity gas. The first 
stars are expected to form in minihalos, M-^^\^ « 10® M0, above the cosmological Jeans 
mass collapsing at 2 ~ 20 — 50 from the highest peaks of the primordial den sity field, 
where cooling is possible by means of molecular hydrogen (Tvir > 10^ K; Tegmark 



et al 19971. Atomic hydrogen cooling has to await the condensation of larger halos. 



with Tvir ^ 10 K and mass ~ 10 M©. Above Mmin the H2 cooling time is shorter 
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than the Hubble time at viriahzation, 
self-gravitating, and stars can form. 



the gas in the central halo regions becomes 



The study of fragmentation of gas in primordial conditions ( Palla et al 2002 1 is 



characterized by some physical simplifications, as compared to present-day star for- 
mation: in the absence of metals and dust the only coolant is hydrogen (atomic and 
molecular), magnetic fields were likely to be dynamically negligible, and the chemistry 
and heating of the primordial gas was not complicated by the presence of a radiation 
background (excepting the cosmic microwave background). 

Simulations of the collapse of primordial molecular clouds ( [Bromm et al[|1999l|2002[ 



Abel et al 2000 Yoshida et al 2006 Gao et al 2007 1 suggest that the first generation 



of stars contained many 'very massive stars' (VMSs) with m*^ > 100 Mq ( Carr et al 
19841. This is because of the slow subsonic contraction - a regime set up by the main 
gas coolant, molecular hydrogen - further fragmentation into sub-components is not 
seen (although it is not clear if this is a numerical effect, rather than due to the gas 



physics, see Glover et al 2008 1 



Moreover, the different conditions of temperature and density of the collapsing 
cloud result in a mass accretion rate over the hydrostatic protostellar core ~ 10^ times 
larger than what observed in local forming stars, suggesting that Pop III stars were 



indeed very massive ( Omukai and Nishi 1998 Ripamonti et al 2002 Tan and McKee 
[2004| ). If the first stars retain their high mass until death, they will collapse after a 
short (~ Myrs) life-time. The final fate depends on the exact mass of the star. Figure 
?? illustrates the fate of primordial stars as a function of their initial mass. 

Between 25 and 140 M0, low-metallicity stars are predicted to form black holes 
directly. The mass of the remnant is about half the star's mass (~ 10 M0 in a 25 M© 
star and ~ 40 Mq in a 100 Mq star; Zhang et al 2008 1. It is likely, however, that if the 
forming black hole is too light, it will not be dynamically stable within the center of 
its host, once stars populate a galaxy. A light black hole might wander within its host, 
dynamically interacting with stars of similar mass, rather than settling at the center 
of the galaxy's potential well. 

Between approximately 140 and 260 Mq lies the domain of pair instability su- 
pernovae. After central helium burning, stars have high enough central entropy that 
they enter a temperature and density regime in which electron/positron pairs are cre- 
ated in abundance, converting internal gas energy into rest mass of the pairs without 
contributing much to the pressure (Barkat et al 1967 Bond et al 1984 1. When this 



instability is encountered, the star contracts rapidly until implosive oxygen and silicon 
burning, depending on the mass of the star, produce enough energy to reverse the 
collapse. These objects are completely disrupted by nuclear-powered explosions. The 
stellar core implodes to a certain maximum temperature that depends on its mass. 



burns fuel explosively, and explodes, leaving no remnants (Kudritzki and Puis 2000 



Fryer et al 2001 1 



In still more massive stars (over 260 Mq on the main sequence), a new phenomenon 
occurs. A sufficiently large fraction of the center of the star becomes so hot that the 
photodisintegration instability is encountered before explosive burning reverses the 
implosion ( [Bond et~al 1984 1 . This uses up all the energy released by previous burning 
stages and, instead of producing an explosion, accelerates the collapse. The nuclear 
energy released by pairs when the star collapses is not sufficient to reverse the implosion 
before the onset of the photodisintegration instability, and the star becomes a black hole 
(Bond et al 1984 Woosley and Weaver 1986 1, sweeping all heavy-element production 
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inside. A massive black hole, containing at least half of the initial stellar mass, is born 
inside the star (Fryer et al 2001 1. 

If VMSs form above 260 Mq, after ~ 2 Myr they would therefore coUapse to MBHs 
with masses intermediate between those of the stellar and supermassive variety. It has 
been suggested by |Madau and Rees| ( |200l] l that a numerous population of MBHs may 
have been the endproduct of the first episode of pregalactic star formation; since they 
form in density high peaks, relic MBHs with mass > 150 M0 would be predicted to 
cluster in the cores of more massive halos formed by subsequent mergers. 

Although this path to MBH formation seems very natural, large uncertainties exist 
on the final mass of PopIII stars. Despite the large number of studies performed, the 
initial mass function of PopIII is still not strongly constrained. As mentioned, it is 
not clear if the very first stars formed in singles or multiples per halo, with the latter 



option leading to a less top-heavy mass function ( Clark et al 2008a Stacy et al 2009 



Turk et al 20091. And, if single objects are formed, simulations of the initial phase of 
formation of these objects only show that massive ~ 10"^ Mq clumps of gas can collapse 
leading to the formation of a very dense, optically thick core of « 0.01 Mq. Gas in the 
envelope accretes into the core, increasing the mass of the protostar. The characteristic 
final mass at the end of the accretion process can be much less than the initial clump 
mass as feedback effects can act, strongly reducing the ability of the core to accrete 
material (McKee and Tan 20081. In addition, the characteristic mass of the star can 
depend on additional external factors as the presence of external UV radiation and the 
temperature of the external cosmic microwave background floor (see the discussion in 
[Trenti et al[ |2009[ ). We do not know if PopIII stars are indeed very massive, and in 
particular if they are above the threshold (~ 260 Mq) for MBH formation. 

It has also been proposed that dark matter, in the form of weakly interacting 
massive particles (WIMPs, e.g. neutralinos) can influence the formation of the first 

If dark matter halos have steep density 



stars (Ripamonti et al 2007 



locco 



20081. 



Navarro et al 



19971, during the formation of the first stars 



profiles (e.g., NFW profile 
the baryonic infall compresses the dark matter further (e.g., adiabatic contraction). The 
central dark matter densities can be high enough that WIMP annihilation can provide 
an extra heat source. The energy released by the annihilation of WIMPs in the stellar 
core overcomes cooling processes and halts stellar collapse. The resulting object, a 'Dark 
Star', is supported by dark matter annihilation rather than nuclear fusion ( Spolyar et al[ 
2008 Freese et al 20081. Dark stars are predicted to be ~ 500 
luminosity (10" 



1000 Mq, with large 



10 Lq), and low surface temperature (<10,000 K). Such massive 
systems would collapse into MBHs, providing alternative MBH seeds. 



2.3 Gas-dynamical processes 



Another family of models for MBH formation relies on the collapse of supermassive 
objects formed directly out of dense gas ( Haehnelt and Rees) 



1993 



Loeb and Rasio 



1994; Eisenstein and Loeb', TM5': 'Bromm and Loeb| |2003r|Koushiappas et a\\ p004 



Begelman et al, ,2006; jLodato and Natarajan 20061. The physical conditions (density, 
gas content) in the inner regions of mainly gaseous proto- galaxies make these locii 
natural candidates, because the very first proto-galaxies were by definition metal-free, 
or at the very least very metal-poor. Enriched halos have a more efficient cooling, 
which in turn favors fragmentation and star formation over the efficient collection of gas 
conducive to MBH formation. It has been suggested that efficient gas collapse probably 
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occurs only in massive halos with virial temperatures Tvir ^ 10 K under metal-free 
conditions where the formation of H2 is inhibited (Bromm and Loeb{ 20031, or for gas 



enriched below the critical metallicity threshold for fragmentation ( Santoro and ShuU 



2006 1. For these systems the tenuous gas cools down by atomic hydrogen only until 
it reaches Tgas ~ 4000 K. At this point the cooling function of the atomic hydrogen 
drops by a few orders of magnitude, and contraction proceeds nearly adiabatically. 

Suppressing molecular-hydrogen formation in metal-free galaxies requires, broadly 
speaking, the presence of dissociating UV radiation. The critical UV fluxes required to 
suppress H2 in massive halos (Tvir > 10*K) are indeed high compared to the expected 
level of the cosmic UV background at the redshifts of interest. [Dijkstra et al| ( |2008[ ) 
suggest that halos prone to H2 suppression are the small subset of all halos (estimated at 
~ 10~®) that sample the bright-end tail of the fluctuating cosmic UV background, due 
to the presence of a close luminous neighbour. Spaans and Silk ( 2006 1 suggest instead 
an alternative route: for isothermally collapsing gas at Tvir ^ 10 K line trapping of 
Lyman a photons causes the equation of state to stiffen, making fragmentation harder, 
in systems with metallicity below ~ 10~^ solar. H2 is naturally destroyed in these 
systems, because of the high gas temperature caused by Lyman a trapping. Finally, it 
has been recently suggested that highly turbulent systems are also likely to experience 
a limited amount of fragmentation, suggesting that efficient gas collapse could proceed 
also in metal-enriched galaxies at later cosmic epochs ( Begelman and Shlosman 2009 1 . 

In such halos where fragmentation is suppressed, and cooling proceeds gradually, 
the gaseous component can cool and contract until rotational support halts the col- 
lapse. In the most common situations, rotational support can halt the collapse before 
densities required for MBH formation are reached. Halos, and their baryonic cores, pos- 
sess angular momentum, J, believed to be acquired by tidal torques due to interactions 
with neighboring halos. This can be quantified through the so-called spin parameter, 

which represents the degree of rotational support available in a gravitational system: 

1/9 5 /2 

Agpin = J\E\ ' /GMj_^ , where E and M/j are the total energy and mass of the halo. 

Let /gas be the gas fraction of a proto-galaxy mass, and fd the fraction of the gas 
that can cool. A mass M = fj^ /gas Mf^ would then settle into a rotationally supported 

2002 1 with a scale radius ~ Aspi„ 



disc ( Mo et al 



1998 



Oh and Haiman 



where 



is the virial radius of the proto-galaxy. Spin parameters found in numerical simulations 
are distributed log-normally in Agpin, with mean Aspin = 0.04 and standard deviation 



a\ — 0.5 (Bullock et al 2001 van den Bosch et al 20021. In a typical high-redshift 
galaxy withT~7">To^ {M^ « lO'' M©; Tvir « 500 pc) the tidally induced angular 
momentum would therefore be enough to provide centrifugal support at a distance 
~ 20 pc from the center, and halt collapse, ultimately leading to the formation of a 
disk. Additional mechanisms inducing transport of angular momentum are needed to 
further condense the gas until conditions fostering MBH formation are achieved. 

Eisenstein and Loeb ( 1995| l and Koushiappas et al ( 2004 \ investigated the formation 



of black holes from low angular momentum material, either in halos with extremely low 
angular momentum, or by considering only the low angular momentum tail of material 
in halos with efficient gas cooling. But even in these models, as in all the others, 
substantial angular momentum transport is required in order for the gas to form a 
central massive object, which ultimately collapses as a result of the post-Newtonian 
gravitational instability. 

An appealing route to efficient angular momentum shedding is by global dynamical 
instabilities, such as the "bars-within-bars" mechanism, that relies on global gravita- 
tional instability and dynamical infall (Shlosman et al 1989 Begelman et al 20061. 
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Self-gravitating gas clouds become bar-unstable when the level of rotational support 
surpasses a certain threshold. A bar can transport angular momentum outward on a 
dynamical timescale via gravitational and hydrodynamical torques, allowing the radius 
to shrink. Provided that the gas is able to cool, this shrinkage leads to even greater 
instability, on shorter timescales, and the process cascades. This mechanism is a very 
attractive candidate for collecting gas in the centers of halos, because it works on a 
dynamical time and can operate over many decades of radius. 

It has been also proposed that gas accumulation in the central regions of protogalax- 
ies can be described by local, rather than global, instabilities. During the assembly of a 
galaxy disc, the disc can become self-gravitating. As soon as the disc becomes massive 
enough to be marginally stable, it will develop structures that will redistribute angular 
momentum and mass through the disc, preventing the surface density from becoming 
too large and the disc from becoming too unstable. To evaluate the stability of the 
disc, the Toomre stability parameter formalism can be used. The Toomre parameter 
is defined as Q = ^gs ' "*^here E is the surface mass density, Cs is the sound speed, 
K = y/2V\^/ R is the epicyclic frequency, and V is the circular velocity of the disc. When 
Q approaches a critical value, Qc, of order unity the disc is subject to gravitational 
instabilities. If the destabilization of the system is not too violent, instabilities lead to 
mass infall instead of fragmentation into bound clumps and global star formation in the 



entire disk (Lodato and Natarajan 20061. This is the case if the inflow rate is below a 

critical threshold Mmax = 2qc^ the disk is able to sustain (where ac ~ 0.12 describes 
the viscosity) and molecular and metal cooling are not important. Such an unstable 
disc develops non-axisymmetric spiral structures, which effectively redistribute angular 
momentum, causing mass inflow (Figure [2] upper panel). This process stops when the 
amount of mass transported to the center is enough to make the disc marginally stable. 
The mass that has to be accumulated in the center to make the disc stable. Ma, is 
obtained by requiring that Q = Qc- This condition can be computed from the Toomre 
stability criterion and from the disc properties, determined from the dark matter halo 



via rvir oc 



2/3 



and angular momentum, via the spin parameter, Agpin, defined 



above (Mo et al 19981 



Ma = /dMh, 



alo 



' 8A 
/dQc 



)( 



J_gas 



1/2 



(2) 



for A < Amax = /dQc/8(/d/jd)(rvir/Tgas)^/^. Here Amax is the maximum halo spin 
parameter for which the disc is gravitationally unstable, is the gas fraction that 
participates in the infall, and is the fraction of the halo angular momentum retained 
by the collapsing gas. Given the mass and spin parameter of a halo, the mass that 
accretes to the center in order to make the disc stable is an upper limit to the mass that 
can go into MBH formation. The efficiency of the mass assembly process ceases at large 
halo masses, where the mass-accretion rate from the halo is above the critical threshold 
for fragmentation and the disc undergoes global star formation instead (Figure|2] lower 
panel). 



Wise et al ( 2008 1 and Regan and Haehnelt ( 2009 1 have performed some idealized 



numerical simulations of the evolution of gaseous discs in haloes with suitable charac- 
teristics. In their simulations, the gas can then lose up to 90% of its angular momentum 
due to supersonic turbulent motions, settling eventually into rotationally supported fat 
discs, in a very compact configuration. These fat, self-gravitating discs are marginally 
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Fig. 2 Color-coded density maps for two exponential gaseous discs, each embedded in a 
dark matter halo, soon after the discs reach the instability regime. Upper panel: in this case, 
fragmentation is suppressed, spiral waves develop in the disc, funnelling gas in the central few 
parsecs. Lower panel: strong cooling allows fragmentation in the disc to set in. 
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stable against gravitational instabilities and resemble the dense structure predicted by 
the analytical calculations. Additional numerical experiments, including all the relevant 
physics, are needed to further explore these possible MBH nurseries. 

After it has efficiently accumulated in the center, the gas made available in the 
central compact region can form a central massive object. The typical masses of gas 
collected within the central few parsecs are of order 10^ — 10^ Mq. We expect that, de- 
pending on how fast and efficiently the mass accumulation proceeds, the exact outcome 
would differ. The gas might form a supermassive star (SMS, mass above ~ 5 x 10^ M©), 
which would eventually collapse and form a black hole. The evolution of a SMS de- 
pends on whether nuclear reactions are taken into account and whether the star has a 
fixed mass or grows via accretion during evolution. It is generally thought that SMSs of 



fixed mass, supported by radiation pressure, will evolve as an n = 3 polytrope (Hoyle 



and Fowler 1963 Baumgarte and Shapiro 1999 Saijo et al''2002^). Unifo rmly rotating 



SMSs cool and contract until they reach the instability poinl[^ Baumgarte and Shapiro 



( 1999 1 investigated the structure and stability of a rapidly rotating SMS in equilibrium. 
As long as the viscous or magnetic braking timescale for angular momentum transfer 



is shorter than the evolutionary timescale of a typical SMS (Zeldovich and Novikov 



|1971[ ), the star will rotate as a rigid body, cool and contract until it reaches the mass 



shedding (where the equatorial angular velocity equals the Kepler velocity). Saijo et al 



( 2002 \ investigated the collapse of a rotating SMS in a post-Newtonian approximation 
to investigate whether the large amount of rotation, causes a rotating SMS to form 
a disk, or if it can collapses directly to a black hole. Rotation seems unable to halt 
the collapse and a MBH is likely to be formed. However, a post-Newtonian calculation 
cannot follow the collapse into the strong-field regime and cannot rigorously address 
the outcome of the collapse. The fate of a marginally unstable, maximally rotating 
SMS of arbitrary mass has been investigated numerically in full general relativity by 



Shibata and Shapiro (20021. They found that the final object is a Kerr-like black hole 



(spin parameter ~0.75) containing 90% of the stellar mass. The fate of an isolated SMS 
is therefore the formation of a MBH. 

If the mass accumulation is fast, however, the outer layers of SMSs are not ther- 
mally relaxed during much of the main sequence lifetime of the star (Begelman 2009 1. 
Gas falls onto the star with increasing specific entropy as a function of time leading 
to entropy stratification. Such SMSs are not necessarily well-represented by n = 3 
polytropes, but rather can have a more complex structure with a convective (poly- 
tropic) core surrounded by a convectively stable envelope that contains most of the 
mass. Hydrogen burning in the core starts when the mass and entropy of the star are 
both relatively low, and adjusts to sustain the star through its more massive stages. 
After exhausting its hydrogen, the core of a SMS will contract and heat up until it 
suffers catastrophic neutrino losses and collapses. The initial black hole, with mass of a 
few M0 , formed as a result of core-collapse subsequently grows via accretion from the 
bloated envelopes that result ('quasistars', an initially low-mass black hole embedded 
within a massive, radiation-pressure-supported envelope; see also [Begelman et al[ 12006) 



2008 1. Over time, the black hole grows at the expense of the envelope, until finally the 



growing luminosity succeeds in unbinding the envelope and the seed MBH is unveiled. 
The key feature of this scenario is that while the black hole is embedded within the 
envelope, its growth is limited by the Eddington limit for the whole quasi-star (this can 



^ We refer the reader to Shapiro 
a SMS of fixed mass to collapse. 



( 2004 I for a comprehensive description of the evolution of 
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be understood in analogy with normal stars, where the central core produces enough 
energy to sustain the whole star, rather than the core only), rather than that appro- 
priate for the black hole mass itself. Very rapid growth can then occur at early times, 
when the envelope mass greatly exceeds the black hole mass. 

The rate at which mass is supplied to the black hole's sphere of influence, Msup, 
can be taken to be proportional to the Bondi rate evaluated at the black hole radius of 
influence. The rate at which matter actually reaches the black hole is suppressed due 



to the back reaction of the energy flux inside the radius of influence ( 


Gruzinov 


1998 


Blandford and Begelman 1999 Narayan et al 2000 Quataert and Gruzinov 


2000 


Mn 



the absence of a wind that removes energy a nd/or angula r momen tum, the accretion 



rate is reduced to Mbh ~ e ^{cs/cfMsnp ([Blandford and Bcgclman| |1999[ [20041 



Comparing the accretion rate to the inflow rate onto the quasistar, the feedback energy 
equals the binding energy of the quasistar before the black hole mass has doubled. Thus 
Bondi accretion, even modified by feedback and a finite rate of angular momentum 
transport, should quickly bring the quasistar to the point where its evolution is driven 
by feedback from the black hole. 

The feedback fiux does not blow apart the quasistar, since this would stop the 
growth of the black hole and therefore the feedback. Instead the quasistar expands 
gradually, allowing the black hole accretion rate to adjust so that the feedback energy 
flux equals the Eddington limit for the instantaneous quasistar mass, A/gfj ~ 2 x 
10~^(£/0.1)~^(M*/10'' Mq)Mq yr"^ The feedback energy flux exceeds the Eddington 
limit for the black hole by a factor M^/MgHi thus, the black hole grows at a super- 
Eddington rate as long as A/* > A/"bh- If the quasistar mass continues to increase at a 
constant rate: 

A/bh ^ lO'^e M0yr"^ (3) 

106 M© ^ ' 

i.e., AfBH oc Af* . As the growth of black hole and quasistar proceed, the photospheric 
temperature of the star decreases, until it reaches a floor minimum temperature (~ 4000 
K) below which no hydrostatic solution for the convective envelope exists, causing the 
convective zone to release radiation at a super-Eddington rate. After a quasistar reaches 
this limiting temperature, it is rapidly dispersed, setting the final limit of the black 
hole seed 

6x 10V4000i^\^^V A'/* \^/^ 



el/2 Tfioor y VlO^^ 



i0 



leaving behind a seed of ~ lO'^M©. The masses of the seeds predicted by different 
models of gas infall and SMS structure vary, but they are typically in the range Mbh ~ 
lO'^ - 10^ Af0 (Figure |4|. 



2.4 Stellar-dynamical processes 



Efficient gas collapse, leading to MBH seed formation, is mutually exclusive with star 
formation, as competition for the gas supply limits the mass available. If star formation 
proceeds in small mini-halos (Tvir < lO'^K), triggered by H2 cooling (Bromm et al 



1999 2002 Abel et al 2000), then by the time more massive halos are built-up, they 



will have been enriched with some traces of metals, brought in by their progenitors. 
Fragmentation and formation of low mass stars starts as soon as gas is polluted by 
metals created in the first PopIII stars (Omukai et al 20081. Formation of 'normal' 
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stars opens up a new scenario for MBH seed formation, if stellar-dynamical rather than 
gas-dynamical processes are at play. This first episode of efficient star formation can 



foster the formation of very compact nuclear star clusters ( Schneider et al 2006 Clark 



et al 2008b I where star collisions can lead to the formation o f a VMS, possibly leaving 

10^ - 10"* Mq (jPevecchi and Volonteri 



a MBH remnant with mass in the range 
[20091 ). 

The possibility that an MBH could in principle form as a result of dynamical 



interactions in dense stellar systems is a long standing idea (|Begelman and Rees 


1978 


Ebisuzaki et al| 2001 


Miller and Hamilton 


2002||Portegies Zwart and McMillan 


2002 


Portegies Zwart et al 


2004 Freitag et al 


2006b|a Giirkanetal 2004 20061. During 



their lifetime collisional stellar systems evolve as a result of dynamical interactions. In 
an equal mass system the central cluster core initially contracts as the system attempts 
to reach a state of thermal equilibrium: energy conservation leads to a decrease in the 
core radius as evaporation of the less bound stars proceeds. As a result the central 
density increases and the central relaxation time decreases. The core then decouples 
thermally from the outer region of the cluster. Core collapse speeds up as it is driven 
by energy transfer from the central denser region ( |Spitzer[ |1987[ ) . This phenomenon 
is greatly enhanced in multi-mass systems like realistic star clusters. In this case the 
gravothermal collapse happens on a shorter timescale as dynamical friction causes the 
more massive stars of mass m, to segregate in the centre on a time-scale ~ /m) 
(where tj^ is the relaxation timescale at the half mass radius, Rjj, and (m) is the mean 
mass of a star in the cluster). If mass segregation sets in before the more massive stars 
evolve out of the main sequence (~ 3 Myr), then a sub-system decoupled from the rest 
of the cluster can form, where star-star collisions can take place in a runaway fashion 



that ultimately leads to the growth of a VMS ( Portegies Zwart et al 1999 1 over a short 
timescale: 



5 X 105 Mq 



(m) 



Gaburov et al ( 2009 1 investigate the process of runaway collisions directly with 
hydrodynamical simulations and show that during collisions large mass losses are likely, 
especially at metallicities larger than lO"'^ solar. Metal-enriched VMSs are therefore 
expected to end their lives as objects less massive than ~ ISOM©, collapsing into MBHs 
with mass ~ few Mq or exploding as pair-instability supernovae. 

The growth of a VMS should be much more efficient at low metallicity. Low metal- 
licity can modify the picture in different ways. First, at sub-solar (but not primordial) 
metallicity, all stars with masses > 4OM0 are thought to collapse directly into a black 



hole without exploding as supernovae (Heger et al 20031. Second, the mass loss due 



to winds is much more reduced in metal-poor stars, which greatly helps in increasing 
the mass of the final remnant. 



Devecchi and Volonteri (20091 investigate the formation of MBHs, remnants of 
VMS formed via stellar collisions in the very first stellar clusters at early cosmic times. 
The main features of their model are as follows: they considered halos with virial 
temperatures Tyir^^lG^ K after the first episode of star formation, hence with a low, 
but non-zero, metallicity gas content. This set of assumptions ensures that (i) atomic 
hydrogen cooling can contribute to the gas cooling process, (ii) a UV field has been 
created by the first stars, and (iii) the gas inside the halo has been mildly polluted by 
the first metals. The second condition implies that at low density H2 is dissociated and 
does not contribute to cooling. The third condition ensures that gas can fragment and 
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Fig. 3 Mean cluster masses (upper left panel), radii (upper right), core collapse timescales 
(lower left) and BH masses (lower right) as a function of redshift (from Devecchi & Volonteri 
2008). The orizontal line in the lower left panel marks the critical core collapse timescale for 
VMS formation. The dashed curves mark the l-cr scatter. 



form low-mass stars only if the gas density is above a certain threshold, ncrit,z (which 
depends on the metalhcity, see |Santoro and Shull[ [2006[ ) : this causes only the highest 
density regions of a proto-galactic disc (see section 2.3) to be prone to star formation. 

In Toomre-unstable proto-galactic discs, such as those described in the previous 
section, instabilities lead to mass infall instead of fragmentation into bound clumps and 
global star formation in the entire disk. The gas inflow increases the central density, and 
within a certain, compact, region n > n^rit.z- Here star formation ensues and a dense 
star cluster is formed. At metallicities ~ 10^'* — 10^^ solar, the typical star cluster 
masses are of order 10^ M0 and the typical half mass radii ~ 1 pc. Most star cluster 
therefore go into core collapse in < 3 Myr, and runaway collisions of stars form a VMS 
(see Section 2.3 and Shapiro 20041, leading to a MBH remnant. Figure |3] summarizes 



the main characteristics of the cluster and MBH populations. 
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As the metallicity of the Universe increases, the critical density for fragmentation 

decreases and stars start to form in the entire protogalactic disk. As gas is consumed 

by star formation, the inflow of gas is no longer efficient, more extended clusters form, 

and the core collapse timescale increases (cfr. Equation [5|. [Devecchi and Volonteri 

(20091 find that typically a fraction ~ 0.05 of protogalaxies at 2; ~ 10 — 20 form black 

hole seeds, with masses 
-3 



1000 — 2000 Mq, leading to a mass density in see d s of a 



Devecchi and Volonteri 



(20091 



few ~ 10 M0/Mpc~ . Most of the assumptions in 
have been conservative, but still the population of seeds is comparable to the case of 
Population III star remnants discussed, for instance, in Volonteri et al (20031. The 



fraction of high-redshift galaxies seeded with a MBH is about a factor of 10 below the 
direct collapse case presented in Volonteri et al (2008b I, where a seed was assumed to 
form with a 100% efficiency whenever a protogalaxy disk was Toomre unstable. 



2.5 Primordial Black Holes 



Primordial black holes may be formed also in the early universe by many processes 



( Carr 2003 and references therein) . Broadly speaking, if within some region of space 
density fluctuations are large, so that the gravitational force overcomes the pressure, 
the whole region can collapse and form a primordial black hole. In the early Uni- 
verse, generically, primordial black holes are formed with masses that roughly equal 



the mass within the particle horizon at the redshift of their formation ( Zel'Dovich and 



Novikov 1967 Hawking 19711. The masses of primordial black holes formed in the 



above mentioned processes range roughly from the Planck Mass (black holes formed 
at the Pl anck epoch) to (black holes formed at the QCD phase transition) up to 
10^ M0 ([Khlopov et ail|2005t. 



Several physical or astrophysical constraints restrict the mass range where primor- 
dial black holes are allowed. Primordial black holes with an initial mass smaller than 
about 5 X 10^^ g are expected to be already evaporated due to Hawking radiation. 
For masses ~ 10^^ g, there are strong bounds from the observed intensity of the dif- 



fuse gamma ray background (Page and Hawking 19761, limiting their contribution to 
the matter density to less than one part in 10*^. For larger masses, constraints can 



be deduced from microlensing techniques (Alcock et al 2000 Tisserand et al 



and from spectral distortions of the cosmic microwave background ( Ricotti et al 
which limit the mass to below ~ 10^ M©. 



2007 



2008 



3 The early growth of massive black holes 



Accretion is inevitable during the "active" phase of a galactic nucleus. Observations tell 
us that AGN are widespread in both the local and early Universe. All the information 
that we have gathered on the evolution of MBHs is indeed due to studies of AGN, 
as we have to await for the Laser Interferometer Space Antenna (LISA) to be able to 
"observe" quiescent MBHs in the distant Universe. 



The classic argument of Soltan (19821, compares the total mass of black holes 



today with the total radiative output by known quasars, by integration over redshift 
and luminosity of the luminosity function of quasars ( |Yu and Tremaine) |2002[ |Elvis| 
et al 2002 Marconi et al 20041. The total energy density can be converted into the 



total mass density accreted by black holes during the active phase, by assuming a 
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mass-to-energy conversion efficiency, e ( Aller and Richstone 2002 Merloni et al 2004 



Elvis et al 2002 Marconi et al 2004 1 



PQSo{< z) = 



/bol(l 



L'${L',z) 



Ho{l + z)y/n„,{l + z)3 + nA 



-.dL'dz 



(6) 



Mace 



~ fholL^ c , is converted into MBH mass 
with /boi the bolometric correction, and e the 



where the mass accretion rate, 
growth, Mbh = /bol(l - 
energy conversion coefficient. 

The similarity of the total mass in MBHs today and the total mass accreted by 
MBHs implies that the last 2-3 e-folds of the mass is grown via radiatively efficient 
accretion, rather than accumulated through mergers or radiatively inefficient accretion. 
However, most ot the e-folds (corresponding to a relatively small amount of mass, say 
the first 10% of mass) could be gained rapidly, such as by, e.g., radiatively inefficient 
accretion. This argument is particularly important at early times, since at least some 
MBHs must have grown rapidly. 

Luminous quasars have indeed been detected at very high redshift, 2 > 6, when 
the Universe was less than 1 Gyr old. Follow-up observations confirmed that at least 
some of these quasars are powered by supermassive black holes with masses ~ 10^ M0 
(Barth et al 2003 Willott et al 20051. How could the seed massive black holes have 



grown rapidly enough within this short timcspan is still an open question (e.g.,| Haiman[ 
[2004, ,Shapiro,, 2005, Volonteri and"Rees, 2056] [Tanaka and Haiman| |20M| ) . 



The rate at which a seed massive black hole can accrete is limited by both 'exter- 
nal' and 'internal' effects. On the one hand, 'external' conditions determine how much 
gas is available, at least initially, to the MBH. Cosmological hydrodynamic simula- 
tions suggest that gravitational collapse produces dense central gas concentrations in 
protogalaxies (e.g., [Pelupessy et al} |2007[ [Wise et al| [2008| |Greif et al[[2008|. Atomic 



densities have been found to reach n ~ 10* cm (e.g. 



Greif et al 



2008 



Bromm and 



Loeb 2003 Wise et al 2008 1. On spatial scales that are resolved in the simulations, gas 



is sufficiently concentrated to enable rapid accretion onto a seed black hole, after the 
radiative feedback from the MBH progenitor star has ceased (if MBHs are the remnant 
of Population HI stars, see [Johnson and Bromm]|2007[ ). 

Given an ample gas supply, rapid accretion might be inhibited by 'internal' feed- 
back, that is the radiative output produced by the newly born quasar itself. |Pelupessy[ 
et al (20071; Johnson and Bromm (20071 suggest that the growth of the MBH is 



severely limited by thermal feedback. However, existing cosmological simulations mod- 
eling MBH growth do not resolve the spatial scales that are needed to treat explicitly 
the radiative processes may alter accretion. The simulations also do not resolve the 
fine structure of the dense gaseous clouds where the MBHs are embedded. Sub-grid 
physics prescriptions are normally adopted for the accretion rate and the radiative 



feedback (e.g., Pelupessy et al 2007 Johnson and Bromm 20071. Recently, however. 



Milosavljevic et al ( 2009 1 have simulated radiatively efficient accretion from a uni- 
form, high-density metal-free protogalactic cloud onto a low-mass seed black hole. 
They found that a compact HH region forms around the black hole causing accretion 
to be intermittent, due to alternating radiation pressure-driven expulsion and external 
pressure-driven fallback. The average accretion rate is 1% of the Bondi accretion rate 
calculated ignoring the radiation's influence, and 30% of the Eddington-limited rate. 
This numerical experiment brings up the question on how black holes can handle 



sprees, that is if/ when the supply is super-critical (Volonteri and Rees 2005 1 the excess 
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radiation is effectively trapped. At very liigfi accretion rate^ radiation pressure cannot 
prevent tlie accretion rate from being super-critical, while the emergent luminosity is 



still Eddington limited in case of spherical or quasi-spherical configurations ( Begelman 



1979 Begelman and Meier 1982"). In the spherical case, though this issue remains 



unclear, it still seems possible that when the inflow rate is super-critical, the radiative 
efficiency drops so that the hole can accept the material without greatly exceeding 
the Eddington luminosity. The efficiency could be low either because most radiation 
is trapped and advected inward, or because the flow adjusts so that the material can 
plunge in from an orbit with small binding energy ( Abramowicz and Lasota 19801. The 
creation of a radiation-driven outflow, which can possibly stop the infall of material, 
is also a possibility. 

The evidence for an accelerated growth of billion solar masses MBH in z>6 galaxies 
is stronger and stronger, though. Ajello et al ( 2009 1 recently published the list of blazars 
detected in the all sky survey by the Burst Alert Telescope (BAT) onboard the Swift 
satellite. The derived luminosity function implies a space density ~ 0.15 Gpc~^ in the 
redshift bin 2 = 3 — 4 of blazars pointing at us with hard X-ray luminosities exceeding 
2 X 10^ erg s^^. Using the IR-optical shape, that is produced by the accretion disk, 
whose radiation is not beamed ( Ghisellini et al [2009 1 find that in all cases of X-ray 
(beamed) luminosity exceeding 2 x lO'" erg s~ , there is an IR-optical (unbeamed) 
spectrum indicating MBHs with masses exceeding 10^ Mq . Since blazars are selected on 
the basis of their X-ray flux relativistically beamed at us, each of these objects implies 
the existence of other ~ 2r^ misaligned similar sources. The space density of MBHs 
having M > 10^ Mq in the z = 3 - 4 redshift range is ~ 68 {r/15 f Gpc~^. The strong 
cosmological evolution revealed by the BAT blazars implies that beamed, jetted sources 
become dominant at the highest MBH masses (> 10^ Mq) and redshifts (z > 5), when 
compared to the cosmic evolution of radio-quiet sources. The number density of blazars 
estimated from the BAT luminosity function is indeed comparable to the total number 
density of > lO^Af0 MBHs at 2; > 6 possibly hosted in massive halos. The presence 
of these astonishing sources suggests that normal 'feedback' may not be at play at the 
highest redshifts. I expect that high-accretion rate events (distinctively possible during 
the violent early cosmic times) will trigger the formation of coUimated outflows (e.g. 
blazars) that do not cause feedback directly on the host (which is pierced through), 
but will deposit their kinetic energy at large distances, leaving the host unscathed (in a 
different context, see |Vernaleo and Reynolds] |2006] ). This is likely if at large accretion 
rates photon trapping decreases the disk luminosity, while concurrently the presence 
of a jet helps dissipating angular momentum, thus promoting efficient accretion. This 
suggestion may explain why high-redshift MBHs can accrete at very high rates without 
triggering self-regulation mechanisms. 



4 Observational signatures of massive black hole seeds 



Figure |4] shows three mass functions for three different MBH 'seed' models: direct 
collapse (|Lodato and Natarajan 2006 1 , runaway stellar mergers in high-redshift clusters 
( [Devecchi and Volonteri 2009 1 , and Population HI remnants ( [Madau and Rees{ |2001[ ) . 



^ We distinguish here between super-Eddington, or super-critical, accretion rate and super- 
Eddington luminosity. A black hole accreting at super-critical values does not necessarily 
radiate at super-Eddingt on luminosity, depending on the properties of the accretion disc 
l |Abramowic2 et al[|1988^ . 
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Fig. 4 Mass function of seed MBHs for three different formation scenarios: direc t collapse 
JVolontcri ct al 2008^ left), runaway stellar mergers in high-redshift c lusters (Dcvecchi andl 
l^lonterij |2009| center), and Population III remnants ( |Madau and Reesl|2001| right). Note the 
different y-axis scale for the Population III case. 



These are the initial conditions that we would like to probe. In this section I will focus 
on two difTerent, extreme, scenarios: 'light seeds', derived from Population III remnants, 
and 'heavy seeds', derived from direct gas collapse. In the following I delineate some 
possible observational tests that can be performed with current and future instruments. 
MBHs have unique properties that make them perfect "beacons" at early cosmic times. 
MBHs that are active, that is, shining as quasars, are the most luminous sources in 
the Universe, making them lighthouses in the early stages of galaxy assembly, when 
the starlight from galaxies is too dim to be detected by our telescopes. Also, black 
holes are not only seen in the electromagnetic spectrum, but they are also sirens for 
gravitational wave detectors, bringing about a novel way of observing our Universe. 



4.1 Tracing MBHs at the earliest times 



As MBHs increase their mass by a large factor during the quasar epoch (z ~ 3 — 4) 
( Yu and Tremaine) 20021, signatures of the seed formation mechanisms are likely more 
evident at earlier epochs. Figure |5] compares the integrated comoving mass density 
in MBHs to the expectations from Soltan-type arguments, assuming that quasars are 
powered by radiatively efficient flows (for details, see [Marconi et al| ( |2004[ |; [Hopkins] 



et al (2007l). The curves differ only with respect to the MBH formation scenario. 
We either assume that seeds are Population III remnants (red curve), or that seeds 
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are formed via direct collapse (Lodato and Natarajan 2006 blue). We study MBH 



evolution within dark matter halos via a Monte-Carlo algorithm based on the extended 
Press-Schechter formalism. The population of MBHs evolves along with their hosts 



according to a "merger driven scenario" , as described in ( Volonteri et al 2003 Volonteri 



and Natarajan 2009). An accretion episode is assumed to occur as a consequence of 



every major merger (mass ratio larger than 1:10) event. During an accretion episode, 
each MBH accretes an amount of mass, AM, that scales with the A/bh ~ relation 
of its hosts (see Volonteri and Natarajan 20091. Accretion starts after a dynamical 



timescale and lasts until the MBH has accreted AM. 

While during and after the quasar epoch the mass densities in the theoretical 
models differ by less than a factor of 2 ('The night in which all cows are black', quoting 
the German philosopher Hegel), at z > 3 the differences become more pronounced. 
Notice that Soltan-type arguments constrain the accreted mass density (dashed curves 
in Figure [5|, rather than the total mass density (solid curves in Figure [5|, making the 
distinction much harder from the observational point of view. 

The initial conditions tend indeed to be erased very fast if accretion is efficient in 
growing MBHs. Volonteri and Gnedin ( 2009[ ) focus on the ultra-high-redshift evolution 
of two distinct populations of seed MBHs: 'light seeds', derived from Population III 
remnants, and 'heavy seeds', derived from direct gas collapse. In both cases they treat 
accretion in a merger-driven fashion (as described above) . As a baseline model, one can 



assume that accretion proceeds at the Eddington rate (see also Salvaterra et al 20071. 
In a second case, accretion during the active phase is based on the extrapolation of the 
empirical distribution of Eddington ratios, A — log(Z/t,ol/^Edd)) found in [Merloni and| 



Heinz ( 2008 1 . We adopt a fitting function of the Eddington ratio distribution as a func- 



tion of MBH mass and redshift Merloni and Heinz ( 2008 1 . We are here extrapolating 



such a model at much higher redshifts and lower MBH masses than originally intended. 
We caution readers in taking the results of this model face value. The main goal of this 
exercise is to probe possible sensible ranges for the accretion rates on MBHs. 

Figure [6] illustrates how, even starting from very different initial mass functions for 
the MBHs (top two panels) hy z — 6 the mass function at the largest masses is similar, if 
accretion is efficient. This washes out any trace of the initial seed population. |Volonteri| 
and Gnedin ( 2009 1 also find that of the contribution to the quasar luminosity budget 
is dominated by MBHs with mass < 10® M©. Such small, low-luminosity MBHs do not 
contribute to the bright end of the luminosity function of quasars, and are therefore 
difficult to account from simple extrapolations of the luminosity function of quasars. 
These small holes are not hosted in extremely massive galaxies residing in the highest 
density peaks (5 to 6~a peaks), but are instead found in more common, "normal" 

systems, ~ 3 a, peaks. Future generation of space-based telescopes, such as JWST 

and IXO, are likely to detect and constrain the evolution of the population of accreting 
massive black holes at early times (z>10). 



4.2 Gravitational waves 



Detection of gravitational waves from seeds merging at the redshift of formation ( Sesana 



et al 20071 is probably one of the best ways to discriminate among formation mech- 



anisms. LISA in principle is sensible to gravitational waves from binary MBHs with 
masses in the range 10"^ — 10® Mq basically at any redshift of interest. A large fraction 
of coalescences will be directly observable by LISA, and on the basis of the detection 
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Fig. 5 Integrated black hole mass density as a function of redshift. Shaded area: constraints 
from Soltan-type arguments, where we have varied the radiative efficiency from a lower limit 
of 6% (upper envelope of the shaded area), to about 30%. Solid lines: total mass density locked 
into nuclear MBHs. Dashed lines: integrated mass density accreted by MBHs. Mode ls based 
on remnants of Population III stars (lowest curve), models based on direct collapse | |Lodato| 
|and Natarajan[|2006^ . Notice that Soltan-type arguments constrain the accreted mass density, 
rather than the total mass density. 



rate, constraints can be put on the MBH formation process. Different theoretical mod- 
els for the formation of MBH seeds and dynamical evolution of the binaries indeed 
predict merger rates that largely vary one from the other (Figure [7|. 



Sesana et al ( 2007 1 and Arun et al ( 2009 1 analyze merger histories of MBHs along 
the hierarchical build-up of cosmic structures (Volonteri et al', 2003 ) . They consider 
two main scenarios for MBH formation (Madau and Rees[ 2001 ,Bcgclman et al, 2006] ) 



Sesana et al ( 2007 1 find that a decisive diagnostic is provided by the distribution of the 



mass ratios in binary coalescences. Models where seeds start large predict that most of 
the detected events involve equal mass binaries. A fraction of observable coalescences, 
in fact, involve MBHs ai z > 10, when MBHs had no time to accrete much mass yet. 
As most seeds form with similar mass, mergers at early times involve MBH binaries 
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Fig. 6 Mass function of MBHs at different cosmic times. Top to bottom: 18 < z < 20; 
16 < 2 < 18; 14 < ^ < 16; 12 < 2 < 14; 10 < 2 < 12; 8 < 2 < 10; 6 < 2 < 8. Left panels: small 
seeds. Right panels: large seed s. Solid histogram: Eddington accretion; dot-dashed histogram: 
Merloni-Heinz accretion. From [Volonteri and Gnedinj l |2009^ . 



with mass ratio ~ 1. In scenarios based on Population III remnants, z > 10 mergers 
involve MBHs with mass below the LISA threshold. The detectable events happen 
at later times, when MBHs have already experienced a great deal of mass growth 
yielding a mass ratio distribution which is flat or features a broad peak at mass ratios 
~ 0.1 — 0.2. This is because both the probability of halo mergers (because of the steep 
host mass function) and the dynamical friction timescale increase with decreasing halo 
mass ratio. Hence, equal mass mergers that lead to efficient binary formation within 
short timescales (i.e., shorter than the Hubble time) are rare, while in more common 
unequal mass mergers it takes longer than an Hubble time to drag the satellite hole 
to the center. Finally, a further, helpful diagnostic of black hole formation models lies 
in the shape of the mass distribution of detected events with S/N>5. Light binaries 
(m < 10^ Mq) are predicted in a fairly large number in Population HI remnant models, 
but are totally absent in direct collapse models. 
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Table 1 Arun et al ( 2009 I evaluated parameter accuracies including also MBH spin infor- 
mation. The four models are referred to as follows. SE: Population III remnant seeds, high 
spins; SC: Population III remnant seeds, low spins; LE: gas-dynamical collapse, high spins; 
LC: gas-dynamical collapse, high spins. For each merger-tree model (SE, SC, LE and LC) we 
list: the total number N of merger events in LISA'S past light cone in a one-year observation; 
the number of events Af^et detectable with SNR larger than 10 in one year; the number for 
which the error in the luminosity distance Dj^ is 10% or less; the number that is localizable 
within 1 and 10 deg^ (^ideg^ ^^'^ ^lOdog^i respectively); and the same for 1 deg^ in angular 
resolution. Results are shown for two different LISA configurations and associated noise curves: 
a "6-link" configuration, which allows the construction of all three independent TDI channels, 
and a "baseline" configuration of 4 links, producing a single Michelson channel. 



Model 


N 






^10deg2 




SE 


80 


33 (25) 


21 (8.0) 


8.2 (1.5) 


2.2 (0.6) 


SC 


75 


34 (27) 


17 (4.4) 


6.1 (0.4) 


1.3 (0.1) 


LE 


24 


23 (22) 


21 (7.7) 


10 (0.8) 


2.2 (0.1) 


LC 


22 


21 (19) 


14 (4.3) 


6.5 (0.5) 


1.8 (0.04) 



Arun et al ( 2009 1 determine the detectability of events by using a code that in- 



cludes both spin precession and higher harmonics in the gravitational-wave signal, and 
carrying out Monte Carlo simulations to determine the number of events that can be 
detected and accurately localized in these population models. LISA will detect a sig- 
nificant fraction of the mergers which occur in the universe: almost all mergers will be 
detected in heavy-seed scenarios, and nearly half of all mergers in small-seed scenarios. 
This is because heavy seed black holes lead to more massive MBH binaries, so they 
can be seen out to larger redshifts. Additionally, heavy-seed models are likely to pro- 
duce more symmetric similar-mass binaries (which produce larger signal-to-noise ratios 
for fixed total mass). In contrast, light-seed models lead to more asymmetric systems. 
Thus, although light-seed MBH mergers could occur more frequently, a smaller frac- 
tion would be observed by LISA due to their smaller total mass and less symmetric 
mass ratios. For heavy-seed models most mergers detectable with high signal-noise- 
ratio (S/N>10) occur in the redshift range 3<2<8, with a peak around z ~ 5 (Figure 
[t]). In the case of fight seeds, mergers are roughly uniform in z over the range 4<z<10. 
Results are summarized in Table 1. 



4.3 Footprints on the Afgu — a relation 

[Volonteri and Natarajan] ( |2009[ ) explore the establishment and evolution of the em- 
pirical correlation between black hole mass (Mbh) a-nd velocity dispersion (a) with 
redshift for two seeding models: 'light seeds', derived from Population HI remnants, 
and 'heavy seeds', derived from direct gas collapse. Even though the seeds themselves 
do not satisfy the Mbh ~ <^ relation initially, the relation can be established and 
maintained at all times if self-regulating accretion episodes are associated with major 
mergers. The massive end of the Mbh ~ <^ relation is established early, and lower mass 
MBHs migrate onto it as hierarchical merging proceeds. How MBHs migrate toward 
the relation depends critically on the seeding prescription. Light seeds initially lie well 
below the Mbk ~ o" relation, and MBHs can grow via steady accretion episodes unhin- 
dered by self-regulation. In contrast, for the heavy seeding model, MBHs are initially 
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redshift 




Fig. 7 Merger rate of MBHs for two different seed models. Top: 'light seeds' from population 
III remnants. Bottom: 'heavy seeds' from gas-dynamical collapse. Black: all mergers. Red: 
mergers detectable with S/N>10 in LISA'S baseline configuration. Blue: mergers detectable 
with S/N>10 in the 6-link configuration. 
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over-massive compared to the empirical correlation, and the host haloes assemble prior 
to kick-starting the growth of the MBH. The slope and scatter of the relation at the 
low-mass end, however, appear to be a consequence of the seeding mechanism and the 
self-regulation prescription. We find that if MBH seeds are massive, ~ 10^ Mq, the 
low-mass end of the A/bh ~ <^ flattens towards an asymptotic value, creating a charac- 
teristic 'plume'. This 'plume' consists of ungrown seeds, that merely continue to track 
the peak of the seed mass function down to late times. 

Since it is during accretion episodes that MBIfs move onto the A/gfj ~ relation, 
AGN are better tracers of the correlation itself, and worse tracers of the original seeds. 
Additionally, differences between seeding models appear only at the low-mass end. 
That's because MBHs move on to the M^j^ — a relation starting from the most massive 
systems at any time, as a consequence of major mergers being more common at high 
redshift for the most massive, biased, galaxies. The implication of this result is that 
flux limited AGN surveys tend to be biased toward flnding MBHs that preferentially 
fall and anchor the Mbh ~ <^ relation - these are systems where the initial conditions, 
the seed mass, have been erased. One of the key predictions is the existence of a large 
population of low mass 'hidden' MBHs at high redshift which are undetectable by 
flux limited AGN surveys. This population of low mass black holes (A/bh < 10*^ AT©) 
are outliers at all epochs on the A/bh ~ <^ relation, and they retain the strongest 
memory of the initial seed mass function, due to their quiet accretion history. More 
accurate measurements of MBH masses below A/bh ~ 10^ Mq will enable us to use 
the measured z = relation to constrain seeding models at high redshift since cosmic 
evolution does not appear to smear out this imprint of the initial conditions. The 
scatter in the observed A/bh ~ relation might also provide insights into the initial 
seeding mechanism (Figure [8|. 



4.4 Massive black holes in low-mass and dwarf galaxies 

The repercussions of different initial efficiencies for seed formation for the overall evo- 
lution of the MBH population stretch from high-redshift to the local Universe. The for- 
mation of seeds in a ylCDM scenario follows the cosmological bias. As a consequence, 
the progenitors of massive galaxies (or clusters of galaxies) have a higher probability 
of hosting MBH seeds (cfr. [Madau and Ree"s]|2001[ ). In the case of low-bias systems, 
such as isolated low-mass galaxies, very few of the high-z progenitors have the deep 
potential wells needed for gas retention and cooling, a prerequisite for MBH formation. 

The signature of the efficiency of the formation of MBH seeds will consequently be 
stronger in low-mass galaxies. Figure |9] (bottom panel) shows a comparison between 
the observed A/bh ~ relation and the one predicted by different models (shown with 
circles), and in particular, from left to right, two models based on the [Lodato and 



[Natarajan (20061 seed masses with Qc ~ 1.5 and 2 and a third model based on lower- 



mass Population HI star seeds. The upper panel of Figure [9] shows the fraction of 
galaxies that do host a massive black hole for different velocity dispersion bins. This 
shows that the fraction of galaxies with a MBH increases with increasing halo masses 
at 2 = 0. A larger fraction of low mass halos are devoid of central black holes for 
lower seed formation efficiencies. This is one of the key discriminants between different 
formation scenarios. 

These arguments can be applied to galaxies in different environments. Observa- 
tionally, the best determination of the MBH (via AGN proxy) population has been 
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Fig. 8 The Mbh — f relation for MBHs at different redshifts. Top: MBHs evolve from an 
initial population of seeds based on the model by Lodato & Natarajan (2006), with Qc = 2 
(the lack of any initial Mbh — f correlation for seeds is clearly seen seen in the far left corner 
of the z = 4 panels, green points). Note the 'plume' of MBHs at cr < 50kms~^ that clearly 
persists even at z = 2 from the earliest epochs. Bottom: A^bh — o" relation for the case of 
Population III remnant seeds. The lack of an initial Mbh ~ correlation for these se eds is also 
evident here and is shown at the bottom of the 2 = 4 panels (green points) . Prom |Volonteri| 
|and NatarajanlpbOQ^ . 
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performed in Virgo. ( [Decarli et al| |2007 1 analyzed nuclear activity in late type galaxies 
in the Virgo cluster using optical data. They conclude that at galaxy mass Mgg^i > 
j^qIO.5 ]y|-^ ^j^g AGN fraction is unity. As a central black hole is a necessary condition 
for AGN activity, we conclude that the black hole occupation fraction must be unity 
as well. Gallo et al ( 2008 1 conducted instead an X-ray survey where they find that the 
incidence of nuclear X-ray activity increases with the stellar mass of the host galaxy: 
only between 3% and 44% of the galaxies with M < 10^^ Mq harbor an X-ray ac- 
tive supermassive black hole. The fracti on rises to betwe e n 49% and 87% in galaxies 



with stellar mass above M < 10 
and 



'Mr. 



Ferrarese et al 



Cote et al 



( 2006| ), 
10^0 



Wehner and Harris 



( 2006[ ) find that, below M* - 
nuclear star clusters, whose mass scales with M« 
massive black holes detected in brighter galaxies ( Magorrian et al 1998 Marconi and 



( 2006 1 

Virgo spheroidals exhibit 
in the same fashion as those of the 



Hunt 2003 Haring and Rix 20041. The existence of a nuclear star cluster does not 



rule out the presence of a MBH, and nuclear star clusters coexist with MBHs over the 



whole range of MBH masses, such as in NGC 3115 (Kormendy et al 19961. Ferrarese 



et al (20061 however suggest that "bright galaxies often, and perhaps always, contain 



supermassive black holes but not stellar nuclei. As one moves to fainter galaxies, nuclei 
become the dominant feature while MBHs might become less common and perhaps 
disappear entirely at the faint end." Our main interest here is to understand if MBHs 
populate preferentially bright galaxies, and if a mass threshold exists for a galaxy to 
host a MBH. 

The record for the smallest known MBH mass belongs to the dw arf Seyfert 1 galax y 



POX 52. It is believed to contain a MBH of mass Mbh ~ 10 Mq (Barth et al 



20041. 



There are also significant non-detections of MBHs in a few nearby galaxies from stellar- 
dynamical observations, most notably the Local Group Scd-type spiral galaxy M33, in 



which the upper limit to Mbh is just a few thousand solar masses (Gebhardt et al 



2001 Merritt et al 20011. Similarly, in the Loc al Group dwarf elliptical galaxy NGC 
205, Mbh < 3.8 x 10* M© |Valluri et ail|2005| . These results suggest that the MBH 
"occupation fraction" in low-mass galaxies might indeed be significantly below unity, 
but at present it is not possible to carry out measurements of similar sensitivity for 
galaxies much beyond the limits of the Local Group. 

Indeed, simple arguments lead us to believe that MBHs might inhabit also the 
nuclei of dwarf galaxies, such as the satellites of the Milky Way and Andromeda, 
today. Indeed, one of the best diagnostics of 'seed' formation mechanisms would be to 
measure the masses of MBHs in dwarf galaxies. As MBHs grow from lower-mass seeds, 
it is natural to expect that a leftover population of progenitor MBHs should also exist 
in the present universe. As discussed above, the progenitors of massive galaxies have 
a high probability that the central MBH is not "pristine", that is, it has increased 
its mass by accretion, or it has experienced mergers and dynamical interactions. Any 
dependence of Mbh on the initial seed mass is largely erased. Dwarf galaxies undergo 
a quieter merger history, and as a result, at low masses the MBH occupation fraction 
and the distribution of MBH masses still retain some "memory" of the original seed 
mass distribution. The signature of the efficiency of the formation of MBH seeds will 
consequently be stronger in dwarf galaxies ( Volontcri ot al^ ^2008b^ . 



Van Wassenhove et al (20101 find that for the most part MBHs hosted in Milky 



Way satellites retain the original 'seed' mass, thus providing a clear indication of what 
the properties of the seeds were. MBHs generated as 'heavy seeds' have larger masses, 
that would favour their identification, their typical occupation fraction is lower, being 
always below 40% and decreasing to less than a % for 'true' dwarf galaxy sizes. Light, 
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Fig. 9 The A/^h— velocity dispersion (ctc) relation at z = 0. Every circle represents the central 
MBH in a halo of given ac- Observation al data are marked by their quoted errorbars, both in 
(Tc, and in A/bh ( |Trem aine et al||2002|. Left to right panels: Qc = 1.5, Qc = 2, Population 
III star seeds. Top panels: fraction of galaxies at a given velocity dispersion which do host a 
central MBH. From | Volonteri et al| ( |2008b| ). 



Population III remnant, seeds have a higher occupation fraction, but their masses 
have not grown much since formation, making their detection harder. The presence of a 
quiescent MBH can be tested dynamically if the region where the gravitational potential 
of the black hole dominates the gravitational potential of the host can be resolved. This 
region is referred to as the sphere of influence of the black hole. We adopt here the 
definition of the sphere of influence as the region within which the enclosed mass in 
stars equals twice the MBH mass. The radius of the sphere of influence is therefore 
defined as: M{r < Ri^f) = 2 x Mbh- 

The lower panel in Figure [lO] plots stellar velocity dispersion against the radius of 
the sphere of influence, -Rinf , estimated for the eight "classical" dwarf spheroidal (dSph) 
satellites of the Milky Way, for which line-of-sight velocities have been measured for 



up to a few thousand stars per galaxy (Walker et al 2009a I. For these objects we 
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Fig. 10 Top panel: number of stars within a given projected radius in Fornax. Lower curve: 
all Fornax memb er stars for which ve locities are currently available in the published kine- 
matic samples of [Walker et al| l |2009a[ | . Upper curve: all Fornax target candidates, including 
unobserved stars, that are sufficiently bright for velocity measurements with existing 6m - 
10m telescopes. Bottom panel: relationship between velocity dispersion (cr) and radius of the 
sphere of influence of MBHs - defined as the sphere that encompasses 2 X Mgu - for ten halo 
realisations. Triangles: we assume that the MBH si ts on the Mb — O" relations hip. Dots: we 
assume a fixed MBH mass, 10^ Mq. Adapted from [Van Wassenhove et al| l [2010| | . 
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adopt the stellar velocity dispersion measurements of Walker et al ( 2009b I , and then 
adopt an MBH mass from the mass-velocity dispersion relationship (|Tremaine et al 



20021. This case describes expectations likely for Pop III remnant seeds. In order to 



calculate the sphere of influence for the dSphs, we consider the best-fitting mass profiles 
from Walker et al ( 2009b 1. In order to evaluate prospects for detecting kinematic 



signatures from such MBHs in real dSphs, the upper panel of Figure [TO] indicates the 
number of spectroscopic target stars within a given projected radius in Fornax, the 
most luminous dSph satellite of the Milky Way. Curves indicate the cumulative surface 
brightness profiles of 1) all Fornax member stars for which velocities are currently 
available in the published kinematic samples of Walker et al. 2009a , and 2) all Fornax 
target candidates, including unobserved stars, that are sufficiently bright for velocity 
measurements with existing 6m - lOm telescopes. The latter profile represents the 
largest samples that are possible at present; unfortunately, these would include fewer 
than 5 stars within the spheres of influence estimated for the classical dSphs. Thus 
even for the brightest nearby dSphs, the detection of any MBH must await the next 
generation of 20-30m telescopes, which may increase kinematic sample sizes by more 
than an order of magnitude. 

Finally, we consider the spheres of influence due to MBH masses of ~ 10^ M0, 
a mass of order of the upper limits derived for the 'massive seed' scenario. For such 
masses, the expected spheres of influence reach ~ 50 pc for the observed dsphs. For 
Fornax, the expected value of R^^^f ~ 30 pc encloses 10 stars in the existing velocity 
sample, and 25 stars in the list of current target candidates. If all these stars are 
observed, the resulting sample may help diagnose whether Fornax has an MBH of 
mass ~ 10^ M0. 

Pushing these limits further and probing the existence of MBHs in dwarf galaxies 
is observationally challenging. An alternative way of detecting MBHs lurking in dwarf 
galaxies would be via their emission, when accreting surrounding material, either from 
a companion star or gas available as recycled material via mass loss of evolved stars 
(Dotti et al. in preparation). A straightforward evidence for MBHs would in fact be 
the presence of AGNs (Greene and Ho 2007 Ho 2008 and references therein). One 



complication in the interpretation of AGN data is the possible contamination by X-ray 



10^ 



lO^M© 



binaries, that have a luminosity comparable to that expected from a 
hole accreting from its surrounding gaseous environment. 

Outside the local group, we can derive an estimate of expected number of dwarf 
galaxies which can possibly host MBHs through two different methods. First, we can 
rely on theoretical models of MBH formation and evolution, and look for the distribu- 



tion of MBHs in dwarf galaxies. Using the dynamical model of Volonteri et al (2008a I, 
we estimate a number density of MBHs, timbh ~ 0.02-0.1 Mpc~'^. 

Second, we can ground our estimate in recent theoretical works that study the 



population of dwarfs as satellites of the Milky Way ( Reed et al 2005 Diemand et al 



2007 Springel et al 20081. These simulations suggest that the number of satellites per 



halo has the following form: 



\ ^host 



«sat 



(7) 



where Usat and «host are the maximum circular velocity of the satellite and the host 
halo, respectively. According to Diemand et al (20071, A''* — 0.021 and a — —3, while 
Springel et all (|2008|) find A^* = 0.052 and a = -3.15. If we extrapolate the Mbh - cr 
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correlation to MBH masses (10 -lO^M©), and assume an isothermal galaxy, then 
«sat ~ 10-20 km s~^. With this formalism we obtain the number of satellites in the 
interesting mass range per dark matter halo (A^sat), where the mass of the halo is 
uniquely determined by its maximum circular velocity. Here, we use A'^ with appropriate 
subscripts to denote occupation number and n to denote number density. The number 
density of dark matter halos can be easily obtained by integrating the modified Press 
& Schechter function iSheth and Tormen 1999) which provides the mass function of 
halos, dn/dAIfi. Therefore we estimate a number density of satellites (per comoving 
cubic Mpc) as: 



Wsat = 



dn 
dM~h 



Ns^t{Mh)dMh = l-3Mpc- 



(8) 



where the lower limit comes from Diemand et al (20071, and the upper limit from 



Springel et al ( 2008 1 . We now have to correct for the fact that not all dwarf galaxies 
are likely to host an MBH, that is, the MBH "occupation fraction" is well below 
unity. To estimate the fraction of dwarfs that host a central MBH, we can rely on the 



models described in the previous paragraph, based on (Volonteri et al 
in which a fraction /mbh 



2003 



0.01-0.1 of dwarfs host an MBH with mass 



2008a I 
- 10""- 



10 Mq . The final estimate for the number density of dwarfs hosting an MBH is then 
^MBH = /MBH'^sat ~ 0.01-0.3 Mpc"'^, in good agreement with the first estimate. 

The detection of gravitational waves from a central MBH in a dwarf galaxy under- 
going a merger is another possible probe. The Einstein Telescope, a proposed third- 
generation ground-based gravitational-wave detector will be able to probe gravitational 
waves in a frequency range reaching down to ~ 1 Hz ( Freise et al 2009 1 . The frequency 
range determines the typical masses of coalescing binaries that could be detected by 
an interferometer; for example, the frequency of gravitational waves emitted from the 
innermost stable circular orbit of a test particle around a Schwarzschild black hole of 
mass M is ~ 4400 Hz(M0/A/). The Einstein Telescope will therefore probe sources 
with masses of hundreds or a few thousand solar masses, which are out of reach of LISA 
or the current ground-based detectors. Since dwarf galaxies have a very quiet merger 
history, we do not expect many MBH-MBH mergers involving dwarf galaxies at the 
present epoch, or in the low-redshift universe. However, gravitational waves may also 
be generated in dwarf galaxies by mergers between the central MBH and stellar rem- 
nants (stellar mass BHs) in the centre of the dwarf. The core stellar densities in ne arby 



dwarf gala xies are typically very low, e.g., t he estimate for Fornax is ~ 10 pc Ma- 



teo 



( 19981) and for Sagittarius is ~ 10 "^pc Majewski et al|(|2005[|. When we calculate 



the event rate of BH-MBH mergers in dwarf galaxies, we have to further correct for 



the fact that only a small fraction of these tiny satellites do indeed form stars ( Bovill 



and Ricotti 2009 and references therein) . Based on Gnedin and Kravtsov ( 2006 1 , we 
estimate that a fraction /* = 0.1 — 0.2 of dwarfs in the Wgat ~ 10 — 20 km s~ range 
formed stars (which will eventually leave behind stellar mass BHs that can merge 
with the central MBH). The number density of MBHs that can be Einstein Telescope 
sources is therefore n^j- = /* riMBH ~ 0.001-0.06 Mpc""^. Therefore, although it is not 
inconceivable that the Einstein Telescope will detect events from dwarf galaxies, any 
events would be serendipitous (^ 1 per year, Gair et al 20091. Detailed calculations 



are needed to understand/prove the robustness of these expectations. 



30 



5 Conclusions 



We can trace the presence of 'super' MBHs at early cosmic times, as the engines pow- 
ering the luminous quasars that have been detected at high redshifts, corresponding 
to about 1 billion years after the Big Bang (Fan et al 2001b I. The HST, Chandra 



and Spitzer satellites, jointly with 8-m class telescopes and large surveys, have made 
important breakthroughs, and observational cosmology probed capable of putting con- 
straints on MBHs, when shining as quasars, up to high redshift (Fan et al 2001b|a 



2004 Barth et al 2003 Willott et al 2003 Walter et al 2004 1 . Some constraints on 



the global accretion history, even at high redshift, can be already put by comparing 



theoretical models predictions to ultra-deep X-ray surveys (Salvaterra et al 20071. The 



early evolution of MBHs, and most notably, what physical mechanism is responsible 
for their formation are however still unknown. We now do know that MBHs are there, 
but we do not know how they got there. 



In this article I focused on three plausible mechanisms of MBH seed formation. 
Broadly speaking, we can divide them into physically-based categories: 'light seeds', 
forming at very early cosmic times {Mgfj ~ 100 — 600 M©, z ~ 20 — 50), 'heavy seeds'. 



forming later on [Mbh - 10 - 10*' Mq, z ~ 

with masses, and at epochs, in between the two previous cases {Mbh — 10"^ M©, z 



10), and 'intermediate seeds', forming 

3i 



10— 15),. Light seeds forming early have a longer time to grow by accretion and MBH- 
MBH mergers, on the other hand, their accretion rates could be depressed in the shallow 
potential wells of the (mini-)halos, especially in the presence of radiative feedback from 
the (mini-)quasar itself (Milosavljevic et al 2009 and references therein). Heavy seeds 



forming later, in more substantial galaxies, are not likely to suffer from the same 
problems, but have had less time to grow. This would be partly compensated by their 
larger initial masses (Volonteri & Begelman 2010). At the current time observational 
constraints are too weak to favor one model against the others, but future X-ray 
missions, such as IXO, and near infrared facilities such as JWST, will have the technical 
capabilities to detect accreting MBHs at z > 6, giving constraints on the accretion 
properties of MBHs at early times. If the mass of the seeds is below ~ 10^ M0, their 
flux is too weak for single sources to be detected electromagnetically. Seeds of mass 
< 10^ Mq can nevertheless be directly identified during their mergers, by detecting 



their emission of gravitational radiation (Hughes 2002 Berti et al 2006 1. Additionally, 



gravitational waves produced during the inspirals of compact objects into MBHs - 
extreme-mass-ratio inspirals (EMRIs) are expected to provide accurate constraints on 
the population of MBHs in the 1O^M0-1O''M0 range (Gair et al. 2010), which is the 
mass range where we can expect 'memory' of the initial conditions, as detailed in 
section 4.4. The combination of electromagnetic and gravitational wave observations 
in the coming years will improve the currently limited constraints on what route, or 
routes, lead to MBH seed formation. 
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